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Abstract 

Hadronic interactions are discussed within the context of the constituent 
quark model. The "Quark Born Diagram" methodology is outlined, ex- 
tensive applications to meson-meson and meson-baryon interactions are 
discussed, and general features of these interactions are highlighted. The 
second half of this document deals with shortcomings of the quark model 
approach and methods to overcome them. These include relativistic 
kinematics, unitarity, nonlocal potentials, coupled channel effects, and 
the chiral nature of the pion. 

1 Ubiquitous Hadronic Interactions 

A microscopic understanding hadron-hadron scattering remains an elusive goal 
of hadronic physics. This is unfortunate because the interactions of hadrons is 
important from a variety of perspectives. At the hadronic level, it provides vital 
insight into the dynamics of quarks and gluons. It is also relevant to the search 
for the quark gluon plasma, since hadronic interactions can mask putative 
signals for the QCD phase transition. Applications extend beyond hadronic 
physics: for example the search for CP violating phases in the final states of 
D of B decays will require correctly accounting for strong interaction final 
state phases|l[]. Of course, hadron-hadron interactions are directly relevant to 
a longstanding goal of nuclear physics - deriving the nuclear force from QCD. 
This issue is not just of intellectual concern since it is important to be able 
to extend our understanding of internuclear forces to extreme conditions (of 
temperature and pressure) so that a variety of astrophysical and cosmological 
issues may be reliably examined. 

Given the importance of the area, it is not surprising that a number of 
techniques have been developed to address the dynamics of hadron-hadron 
interactions. These roughly fall into two classes: those which treat hadrons as 
elementary fields and those which attempt to describe the interactions using 



QCD as the starting point. Among the former are potential approaches^ 
(Bonn, Paris, Argonne), relativistic hydrodynamic approaches]^, and a variety 
of effective field theories[Q (chiral perturbation theory, effective NN theory). 
Among the latter are lattice gauge theory]^, Schwinger-Dyson models^], light 
cone field theory]^, a multigluon dipole interaction model 0, and constituent 
quark models. 

Unfortunately, it is not possible to summarize the advantages and faults of 
all of these methods here. Rather, the remainder of this document focuses on 
attempts to describe hadronic interactions which are based on the nonrelativis- 
tic constituent quark model (CQM). The CQM is beyond a doubt the most 
widely used description of the static properties of hadrons - largely because it 
is able to describe hundreds of experimental data with a handful of parame- 
ters in a comparatively simple picture. Thus it is no surprise that attempts to 
describe hadronic interactions with the CQM have a long history Despite 
the well-known short comings of the CQM (some of which will be discussed 
below) the benefits are immediately apparent: the CQM may be applied to 
any hadronic interactions and the predictions are essentially parameter-free. 
By this I mean that the parameters of the CQM are very strongly fixed by 
comparison with static hadron properties - there is no wiggle room when com- 
puting interactions! Of course the extension of the CQM to dynamic properties 
of hadrons requires some extrapolation^; however, this should not form a bar- 
rier in itself. To paraphrase Feynman, "Trust your model and see how far it 
takes you." 



2 Quark Born Diagrams 



In the following we shall consider a CQM which includes Coulomb and linear 
central potentials, a spin-spin colour hyperfine interaction, and possibly spin- 
orbit and tensor interactions. It may be somewhat of a surprise that the Breit- 
Fermi interactions are included here since they are normally only required to 
achieve detailed agreement with spectroscopic and other static properties of 
the hadrons. However, as we shall see shortly, subleading (in v/c) interactions 
can dominate hadronic scattering! 

A great deal of effort has been expended on variational |TD[ and resonating 



group approaches [|r^] to scattering in the CQM. Here I describe a simplified 
approach where one evaluates the T- matrix at Born order [|12|, |13|, called the 



"Quark Born Diagram" (QBD) method. The T-matrix for meson-meson scat- 

"'^For example, the dynamics of flux tubes may be essentially ignored in conventional 
meson and baryons. This is no longer the case for multiquark systems. 



tering is shown diagrammatically in Fig 1 . The diagrams represent momentum 
flow in the Born order term of the Neumann series and must be attached to ex- 
ternal mesonic wavefunctions. Quark exchange must occur to maintain asymp- 
totic colour singlet states. In doing these computations one needs to be aware 
that scattering with composite objects can be very different than pointlike ob- 
jects. For example, hermiticity of the scattering amplitude is no longer trivial 
but requires that the wavefunctions be exact eigenstates of Hq. This property 
is related to the "post-prior" discrepancy ||T3|, M. Furthermore, consistency of 



the kinematical relationships must be maintained with the Hamiltonian (non- 
relativistic Hamiltonians require nonrelativistic kinematics to maintain her- 
miticity). Also, one cannot change parameters at will! For example, mesonic 
radii or quark masses cannot be independently modified. One must instead ad- 
just a Hamiltonian parameter, solve for the spectrum again, and then evaluate 
the T-matrix. 
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Fig. 1: Diagrams Contributing to Meson-Meson Scattering 
2.1 Applications 

The application of the Quark Born formalism requires some care. For exam- 
ple, if one is to restrict attention to the terms shown in Fig 1, then channels 
where strong quark annihilation effects are expected (such as J = vrvr) should 
be avoided. Furthermore, if the predicted interactions are strong, the Born 
order results should be iterated (this is discussed below). The Quark Born 
Diagram method has been applied to a variety of hadronic reactions. These 
include 1=2 vrvr scattering [|12|, [T5[| (with surprising agreement considering the 
relativistic and chiral nature of the pion), 1=3/2 ttK scattering [jl6[ (testing 
Bose symmetry breaking due to the strange quark mass), KN scattering ||17|| 
(demonstrating surprising agreement in the S-wave and a dramatic failure in 



the P-wave[|T^]), short range A^A^ scattering [|1^] (in agreement with resonating 
group computations), BB scattering pOf (in agreement with lattice computa- 
tions), J/ip — TT scattering [^] (in strong disagreement with previous estimates 
and in agreement with rudimentary data), vrp scattering [|l5l (an examination 
of the generation of hadronic spin orbit and tensor interactions from the quark 
level), and possible meson-meson bounds states||2^ (the /o(1710) may be iden- 
tified K*K* — uj(j) bound state). 

A comparison of the predicted and experimental isotensor S-wave tttt phase 
shifts is shown in Fig 2. A sample cross section prediction relevant to charm 
suppression at RHIC is shown in Fig 3. 
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Fig. 2: The / = 2 S-wave tttt Phase Shift. The sohd hne is the QBD prediction |2|, 
the dashed hne is the one loop chiral prediction while the lower dashed line is 
the tree order chiral result. Data are from IM. 




3 Limitations and Extensions of the CQM 



As mentioned above, the CQM in general, and the QBD formahsm in par- 
ticular, suffers from several inadequacies. We summarize several of these and 
discuss methods for addressing them. 



3.1 Unitarity, Relativistic Kinematics, and Nonlocality 

Perhaps the simplest problem arises when the scattering is so strong that the 
tree level diagrams of the QBD method are inaccurate. This may be tested by 
comparing QBD predictions to more complete resonating group calculations 
(which iterate the scattering to all orders and can, in principle, include the 
effects of coupled channels, wavefunction distortion, etc). Such a comparison 



of pp scattering in isospin 2 was made in Ref. ||13[ where the accuracy of the 
QBD results were explicitly demonstrated. Nevertheless, it is worth observing 
that strong interactions may be accounted for by extracting effective poten- 
tials for the process in question and iterating the potential in the appropriate 
Schrodinger equation. An example of this is shown in Fig 4, where the effective 
TTTT interaction has been extracted from the QBD tttt T-matrix and iterated in 
the nonrelativistic Schrodinger equation. This is shown as a dashed line. Evi- 
dently, the agreement with data is ruined by this procedure (I am presenting 
the worst possible case - the procedure works very well in general). This is 
because the light mass of the pions (and the large invariant masses at which 
the formalism is being applied) requires relativistic kinematics to be employed 
(this has been used in the QBD prediction of Fig 2). Thus the nonrelativistic 
Schrodinger equation is inappropriate for iterating the effective vrvr interaction. 
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Fig. 4: 1=2 TTTT Scattering. The solid line is the QBD prediction, the dashed line is the 
nonrelativistic local prediction, while the dotted line is the nonlocal, relativistic prediction. 



A related problem is the nonlocality of the vrvr T-matrix (because composite 
particles are being scattered). Any extracted effective potential which is local 



will induce unknown errors after being iterated. Fortunately, all three problems 
may be dealt with by choosing to solve for the T-matrix directly in momentum 
space: 

TE{k',k) = V{k\k) + J d'pV{k',p)^-^^-^TE{p,k) (1) 

As is evident from the equation, nonlocal effective potentials may be directly 
employed (in fact these are proportional to the QBD T-matrices), unitarity 
is automatically restored, and nonrelativistic kinematics may be incorporated 
by use of the relativistic dispersion relation for E{p). This is illustrated as the 
dotted line in Fig 4, where one sees that the effect of iterating the tttt potential 
is to weaken the scattering at higher invariant mass. 



3.2 Annihilation to Hybrids 

To this point attention has been restricted to 'exotic' channels such as 1=2, 
where resonance contributions are forbidden. We now examine the issues in- 
volved in relaxing this constraint. Two possible intermediate states may be 
realized (at least at lowest order in the Fock space expansion), annihilation to 
intermediate mesons or annihilation to intermediate hybrid mesons. The latter 
process involves an intermediate state consisting of a quark - antiquark pair 
and a 'gluon' (where the gluon may be an excited flux tube or a constituent 
gluon, depending on one's picture of soft glue). 

The nonrelativistic reduction of the one gluon exchange potential which 



describes the coupling of a qq pair to an intermediate perturbative gluon is|25 



K„„ = ^(2 + S,.S,)^(n4|. (2) 

However, the intermediate state in this case is a hybrid, which typically has 
a mass some 1000 MeV above low lying mesons with the same quark content. 
Thus one may expect that it is more realistic to employ a gluon propagator 
with a fictitious mass of roughly 1 GeV. One may incorporate this into the ex- 
pression above by multiplying it by a factor / which is to be fit to the data and 
which we expect to be roughly -1. This rather speculative adjustment can be 
verified by comparing reactions with no annihilation to similar reactions where 
annihilation is permitted. For example, 1=2 tttt scattering may be compared 
to 1=0 tttt; 1=3/2 Kn scattering may be compared to 1=1/2 Ktt; and K^N 
scattering may be compared to vrA^ scattering. All indicate that a negative 
value of / is required; a fit yields / ~ — 2.6p5|. 



3.3 Coupled Channels 



Even a cursory examination of typical hadronic scattering data reveals the 
importance of intermediate resonance states. Unfortunately the mechanism 
by which hadrons couple is poorly understood and surely involves compli- 
cated nonperturbative gluodynamics. The current best guess is the purely 
phenomenological ^Pq model in which qq pairs are created with vacuum quan- 
tum numbers and combine with the parent quarks to produce the daughter 
mesons. Extensive calculations of meson and baryon decays have been made 
with moderate success (typical errors in the amplitude are 20% or less)||26[]. 

Incorporating the ^Pq model directly into the quark model is the most direct 
way to include the effects of intermediate resonances. This may be achieved 
most simply by writing the quark model in second quantized notation and 
including a ^Pq term which creates and annihilates qq pairs. This is multiplied 
by a constant 7 which may be determined by comparison to a specific channel 
(say, p^tttt). 

H = dx{-- — h%-- — 44)+7 / c?a;(6i,(T- V4 + H.C.). 
J 2mq 2mg J 

+ ^Jdxdy {blbl + dldl)V{x - y){byb^ + dyd^). (3) 

The field theory is simplified by restricting the Fock space to the meson and 
meson-meson sectors of interest. Thus we project onto \BC) by making 
the following Ansatz for the exact eigenstate: 

1^) = J^Ain-r2)b\dl\Q) 

+ / E ^Bc C'^'''''''~''' )Mri - rs)Mr2 - r,)b\4bldl\0) (4) 

BC ^ 

Varying the reduced Hamiltonian matrix element with respect to the unknown 
meson (p^ and meson-meson ^bc wavefunctions yields the coupled channel 
Schrodinger equation: 

EifA{r) = Hgg{r)(pA{r) 
--ijt-iVB + Vc + VBc)^m{r/2-x)^^c{r/2 + x)^Bc{-r/2), (5) 

— + f f KE{x,y,R)^BciR')+ f fvE{x,y,R)<ilBc{R') 



- S ■ (Vb + Vc + Vbc)^qbVocVa{-'2R) 

= E^BciR) +eJ Ne{x, y, R)^Bc{R') (6) 

Here r is the interquark radius in the meson channel and R is the intermeson 
distance in the meson-meson channel. Remarkably there is a simple rela- 
tionship between these coordinates: R = r/2. Ke, Ve, and Ne represent 
the exchange kinetic energy, potential, and normalization kernels respectively. 
Wavef unctions with a 'naught' subscript (fo represent mesonic wavef unctions 
without the effects of channel mixing (so that we have assumed no wavefunc- 
tion distortion in deriving this equation). The first of these equations is the 
nonrelativistic quark model (Hgg) supplemented with a term which couples it 
to the meson-meson continuum - thereby 'unquenching' the quark model. The 
second equation is the resonating group equation which describes meson-meson 
scattering in the CQM (the Born order T-matrix for this equation is provided 
by the QBD). The term proportional to 7 provides the desired coupling to 
intermediate resonances. 

Eqns (5,6) may be solved with standard coupled channel methods and 
the effects of unquenching the quark model and of intermediate resonances in 
scattering problems may be examined. An example of this is given in Fig 5 
where the effect of coupling virtual ps to the vrvr P-wave channel is studied. 
The p mass has been shifted down 80 MeV while the bare width is 10 MeV, 
this corresponds to an 'RPA' width of 90 MeV[0]. 

40 I 1 1 1 1 1 1 1 1 

30 - 
20 - 

10 - - 




-10 - + - 

-20 - 
-30 - 

_40 I 1 1 1 1 1 1 1 

0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

m_pipi (GeV) 

Fig. 5: 1=1 TTTT Scattering. The arrows indicate the locations of the 'bare' p and p' mesons. 
3.4 Chiral Pions 

The apparently successful description of vrvr scattering evident in Fig 1 is per- 
haps surprising given the chiral nature of the pion and its interactions at low 



energy. Explaining this success will go a long way towards explaining the un- 



warranted success of the CQM in the light quark sector |2^. One way to do 
this is to construct a model of strong QCD which incorporates the physics of 
chiral symmetry breaking at a microscopic level. In this way one may compute 
TTTT scattering with composite particles while observing the dictates of chiral 
symmetry. It is possible to construct such a model by assuming a nontrivial 
QCD vacuum (typically a BCS-type vacuum) and building states on this vac- 
uum with the random phase approximation to the full Bethe Salpeter equation. 
This has been done in Coulomb gauge QCD[^, with similar calculations 
in the Schwinger-Dyson approach pl| 



4 Conclusions 

Studying and understanding hadronic interactions is vital to hadronic physics 
and is an important part of nuclear and electroweak physics. It is probable that 
a microscopic description of hadronic interactions is necessary if one wishes 
to understand these phenomena in extreme conditions or in poorly known 
channels. The constituent quark model provides an excellent starting point for 
developing the understanding required to construct a reliable model of strong 
QCD. In the meantime, it also serves as an excellent phenomenological guide to 
the interpretation of scattering experiments. The development of continuum 
field theoretic models capable of describing hadrons and their interactions is in 
its infancy - we look forward to their maturation and application to reaction 
processes. 
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tal. A large portion of the work described here was performed with Ted Barnes 
and Adam Szczepaniak. The author is grateful to Paul Geiger for discussions 
on coupled channels. This work was supported by the US DOE under contracts 
DE-FG02-00ER41135 and DE-AC05-84ER40150. 
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Abstract 

Hadronic interactions are discussed within the context of the constituent 
quark model. The "Quark Born Diagram" methodology is outlined, ex- 
tensive applications to meson-meson and meson-baryon interactions are 
discussed, and general features of these interactions are highlighted. The 
second half of this document deals with shortcomings of the quark model 
approach and methods to overcome them. These include relativistic 
kinematics, unitarity, nonlocal potentials, coupled channel effects, and 
the chiral nature of the pion. 

1 Ubiquitous Hadronic Interactions 

A microscopic understanding hadron-hadron scattering remains an elusive goal 
of hadronic physics. This is unfortunate because the interactions of hadrons is 
important from a variety of perspectives. At the hadronic level, it provides vital 
insight into the dynamics of quarks and gluons. It is also relevant to the search 
for the quark gluon plasma, since hadronic interactions can mask putative 
signals for the QCD phase transition. Applications extend beyond hadronic 
physics: for example the search for CP violating phases in the final states of 
D of B decays w'ill require correctly accounting for strong interaction final 
state phases[l]. Of course, hadron-hadron interactions are directly relevant to 
a longstanding goal of nuclear physics - deriving the nuclear force from QCD. 
This issue is not just of intellectual concern since it is important to be able 
to extend our understanding of internuclear forces to extreme conditions (of 
temperature and pressure) so that a variety of astrophysical and cosmological 
issues may be reliably examined. 

Given the importance of the area, it is not surprising that a number of 
techniques have been developed to address the dynamics of hadron-hadron 
interactions. These roughly fall into two classes: those which treat hadrons as 
elementary fields and those which attempt to describe the interactions using 



QCD as the starting point. Among the former are potential approaches [2] 
(Bonn, Paris, Argonne), relativistic hydrodynamic approaches [3], and a variety 
of effective field theories[4] (chiral perturbation theory, effective NN theory). 
Among the latter are lattice gauge theory[5], Schwinger-Dyson models[6], light 
cone field theory [7], a multigluon dipole interaction model [8], and constituent 
quark models. 

Unfortunately, it is not possible to summarize the advantages and faults of 
all of these methods here. Rather, the remainder of this document focuses on 
attempts to describe hadronic interactions which are based on the nonrelativis- 
tic constituent quark model (CQM). The CQM is beyond a doubt the most 
widely used description of the static properties of hadrons - largely because it 
is able to describe hundreds of experimental data with a handful of parame- 
ters in a comparatively simple picture. Thus it is no surprise that attempts to 
describe hadronic interactions with the CQM have a long history [9]. Despite 
the well-known short comings of the CQM (some of which will be discussed 
below) the benefits are immediately apparent: the CQM may be applied to 
ayiy hadronic interactions and the predictions are essentially parameter-free. 
By this I mean that the parameters of the CQM are very strongly fixed by 
comparison with static hadron properties - there is no wiggle room when com- 
puting interactions! Of course the extension of the CQM to dynamic properties 
of hadrons requires some extrapolation^; however, this should not form a bar- 
rier in itself. To paraphrase Feynman, "Trust your model and see how far it 
takes you." 

2 Quark Born Diagrams 

In the following we shall consider a CQM which includes Coulomb and linear 
central potentials, a spin-spin colour hyperfine interaction, and possibly spin- 
orbit and tensor interactions. It may be somewhat of a surprise that the Breit- 
Fermi interactions are included here since they are normally only required to 
achieve detailed agreement with spectroscopic and other static properties of 
the hadrons. However, as we shall see shortly, subleading (in v/c) interactions 
can dominate hadronic scattering! 

A great deal of effort has been expended on variational [10] and resonating 
group approaches[ll] to scattering in the CQM. Here I describe a simplified 
approach where one evaluates the T-matrix at Born order [12, 13], called the 
"Quark Born Diagram" (QBD) method. The T-matrix for meson-meson scat- 

^For example, the dynamics of flux tubes may be essentially ignored in conventional 
meson and baryons. This is no longer the case for multiquark systems. 



tering is shown diagrammatically in Fig 1. The diagrams represent momentum 
flow in the Born order term of the Neumann series and must be attached to ex- 
ternal mesonic waveftmctions. Quark exchange must occur to maintain asymp- 
totic colour singlet states. In doing these computations one needs to be aware 
that scattering with composite objects can be very difl^erent than pointlike ob- 
jects. For example, hermiticity of the scattering amplitude is no longer trivial 
but requires that the wavefunctions be exact eigenstates oi Hq. This property 
is related to the "post-prior" discrepancy[13, 14]. Furthermore, consistency of 
the kinematical relationships must be maintained with the Hamiltonian (non- 
relativistic Hamiltonians require nonrelativistic kinematics to maintain her- 
miticity). Also, one cannot change parameters at will! For example, mesonic 
radii or quark masses cannot be independently modified. One must instead ad- 
just a Hamiltonian parameter, solve for the spectrum again, and then evaluate 
the T-matrix. 




Fig. 1: Diagrams Contributing to Meson-Meson Scattering 
2.1 Applications 

The application of the Quark Born formalism requires some care. For exam- 
ple, if one is to restrict attention to the terms shown in Fig 1, then channels 
where strong quark annihilation effects are expected (such as 7 = tttt) should 
be avoided. Furthermore, if the predicted interactions are strong, the Born 
order results should be iterated (this is discussed below). The Quark Born 
Diagram method has been applied to a variety of hadronic reactions. These 
include 1=2 tttt scattering[12, 15] (with surprising agreement considering the 
relativistic and chiral nature of the pion), 1=3/2 ttK scattering[16] (testing 
Bose symmetry breaking due to the strange quark mass), KN scattering[17] 
(demonstrating surprising agreement in the S-wave and a dramatic failure in 



the P-wave[18]), short range NN scattering[19] (in agreement with resonating 
group computations), BB scattering [20] (in agreement with lattice computa- 
tions), J/il) — TT scattering[21] (in strong disagreement with previous estimates 
and in agreement with rudimentary data), vrp scattering[15] (an examination 
of the generation of hadronic spin orbit and tensor interactions from the quark 
level), and possible meson-meson bounds states[22] (the /o(1710) may be iden- 
tified as a K*K* — uj(j) bound state). 

A comparison of the predicted and experimental isotensor S-wave tttt phase 
shifts is shown in Fig 2. A sample cross section prediction relevant to charm 
suppression at RHIC is shown in Fig 3. 
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Fig. 2: The 1 = 2 S-wave tttt Phase Shift. The solid Hne is the QBD prediction[12], 
the dashed Hne is the one loop chiral prcdiction[23], while the lower dashed line is 
the tree order chiral result. Data are from [24]. 




3 Limitations and Extensions of the CQM 



As mentioned above, the CQM in general, and the QBD formalism in par- 
ticular, suffers from several inadequacies. We summarize several of these and 
discuss methods for addressing them. 

3.1 Unitarity, Relativistic Kinematics, and Nonlocality 

Perhaps the simplest problem arises when the scattering is so strong that the 
tree level diagrams of the QBD method are inaccurate. This may be tested by 
comparing QBD predictions to more complete resonating group calculations 
(which iterate the scattering to all orders and can, in principle, include the 
effects of coupled channels, wavefunction distortion, etc). Such a comparison 
of pp scattering in isospin 2 was made in Ref. [13] where the accuracy of the 
QBD results were explicitly demonstrated. Nevertheless, it is worth observing 
that strong interactions may be accounted for by extracting effective poten- 
tials for the process in question and iterating the potential in the appropriate 
Schrodinger equation. An example of this is shown in Fig 4, where the effective 
TTTT interaction has been extracted from the QBD tttt T-matrix and iterated in 
the nonrelativistic Schrodinger equation. This is shown as a dashed line. Evi- 
dently, the agreement with data is ruined by this procedure (I am presenting 
the worst possible case - the procedure works very well in general). This is 
because the light mass of the pions (and the large invariant masses at which 
the formalism is being applied) requires relativistic kinematics to be employed 
(this has been used in the QBD prediction of Fig 2). Thus the nonrelativistic 
Schrodinger equation is inappropriate for iterating the effective tttt interaction. 
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Fig. 4: 1=2 TTTT Scattering. The solid line is the QBD prediction, the dashed Hne is the 
nonrelativistic local prediction, while the dotted line is the nonlocal, relativistic prediction. 



A related problem is the nonlocality of the tttt T-matrix (because composite 
particles are being scattered). Any extracted effective potential which is local 



will induce unknown errors after being iterated. Fortunately, all three problems 
may be dealt with by choosing to solve for the T-matrix directly in momentum 
space: 

TE{k',k) = V{k',k) + I d'pV{k',p)^-^^-^TE{p,k) (1) 

As is evident from the equation, nonlocal effective potentials may be directly 
employed (in fact these are proportional to the QBD T-matrices), unitarity 
is automatically restored, and nonrelativistic kinematics may be incorporated 
by use of the relativistic dispersion relation for E{p). This is illustrated as the 
dotted line in Fig 4, where one sees that the effect of iterating the tttt potential 
is to weaken the scattering at higher invariant mass. 

3.2 Annihilation to Hybrids 

To this point attention has been restricted to 'exotic' channels such as 1=2, 
where resonance contributions are forbidden. We now examine the issues in- 
volved in relaxing this constraint. Two possible intermediate states may be 
realized (at least at lowest order in the Fock space expansion), annihilation to 
intermediate mesons or annihilation to intermediate hybrid mesons. The latter 
process involves an intermediate state consisting of a quark - antiquark pair 
and a 'gluon' (where the gluon may be an excited flux tube or a constituent 
gluon, depending on one's picture of soft glue). 

The nonrelativistic reduction of the one gluon exchange potential which 
describes the coupling of a qq pair to an intermediate perturbative gluon is[25] 

K„„ = ^(2 + S<.S,)i(f,,)||. (2) 

However, the intermediate state in this case is a hybrid, which typically has 
a mass some 1000 MeV above low lying mesons with the same quark content. 
Thus one may expect that it is more realistic to employ a giuon propagator 
with a fictitious mass of roughly 1 GeV. One may incorporate this into the ex- 
pression above by multiplying it by a factor / which is to be fit to the data and 
which we expect to be roughly -1. This rather speculative adjustment can be 
verified by comparing reactions with no annihilation to similar reactions where 
annihilation is permitted. For example, 1=2 tttt scattering may be compared 
to 1=0 tttt; 1=3/2 Ktt scattering may be compared to 1=1/2 Kn; and K'^N 
scattering may be compared to 'wN scattering. All indicate that a negative 
value of / is required; a fit yields / ~ —2. 6 [2 5]. 



3.3 Coupled Channels 



Even a cursory examination of typical hadronic scattering data reveals the 
importance of intermediate resonance states. Unfortunately the mechanism 
by which hadrons couple is poorly understood and surely involves compli- 
cated nonperturbative giuodynamics. The current best guess is the purely 
phenomenological ^Pq model in which qq pairs are created with vacuum quan- 
tum numbers and combine with the parent quarks to produce the daughter 
mesons. Extensive calculations of meson and baryon decays have been made 
with moderate success (typical errors in the amplitude are 20% or less) [26]. 

Incorporating the ^Pq model directly into the quark model is the most direct 
way to include the effects of intermediate resonances. This may be achieved 
most simply by writing the quark model in second quantized notation and 
including a ^Pq term which creates and annihilates qq pairs. This is multiplied 
by a constant 7 which may be determined by comparison to a specific channel 
(say, p^-k-k). 



dx ( - - — bU^ - - — dU^) +'y dx (6j.cr • V4 + H.c). 

ZTTlg ZTJlq J 

+ ^ J dxdy{blbl + dldl)V{x - y){byb^ + dyd^). (3) 

The field theory is simplified by restricting the Fock space to the meson and 
meson-meson sectors of interest. Thus we project onto \A), \BC) by making 
the following Ansatz for the exact eigenstate: 

I*) = I ipA{r,-r2)b\dl\0) 

+ ji: ^Bc C'^''~''~'' )Mri - r,)Mr2 - n)bUlbUl\0) (4) 



BC 



Varying the reduced Hamiltonian matrix element with respect to the unknown 
meson cpA and meson-meson wavefunctions yields the coupled channel 
Schrodinger equation: 

E(fA{r) = Hqq{r)(pA{r) 
-7|s-(Vs + Vc + ^Bc)^oB{r/2-x)(poc{r/2 + x)^Bc{-r/2), (5) 

-V^ + 1 1 KE{x,y,R)^Bc{R') + I lvE{x,y,R)^Bc{R') 



2/^13, 
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- 87 y S • (Vb + Vc + V Bc)'PQB'PQC'PA{-'iR) 

= E^Bc{R) +E j Ne{x, y, R)^BciR') (6) 

Here r is the interquark radius in the meson channel and R is the intermeson 
distance in the meson-meson channel. Remarkably there is a simple rela- 
tionship between these coordinates: R = r/2. Ke, Ve, and Ne represent 
the exchange kinetic energy, potential, and normalization kernels respectively. 
Wavefunctions with a 'naught' subscript (po represent mesonic wavefunctions 
without the effects of channel mixing (so that we have assumed no wavefunc- 
tion distortion in deriving this equation). The first of these equations is the 
nonrelativistic quark model {Hgq) supplemented with a term which couples it 
to the meson-meson continuum - thereby 'unquenching' the quark model. The 
second equation is the resonating group equation which describes meson-meson 
scattering in the CQM (the Born order T-matrix for this equation is provided 
by the QBD). The term proportional to 7 provides the desired coupling to 
intermediate resonances. 

Eqns (5,6) may be solved with standard coupled channel methods and 
the effects of unquenching the quark model and of intermediate resonances in 
scattering problems may be examined. An example of this is given in Fig 5 
where the effect of coupling virtual ps to the tttt P-wave channel is studied. 
The p mass has been shifted down 80 MeV while the bare width is 10 MeV, 
this corresponds to an 'RPA' width of 90 MeV[27]. 
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Fig. 5: 1=1 TTTT Scattering. The arrows indicate the locations of the 'bare' p and p' mesons. 
3.4 Chiral Pions 

The apparently successful description of tttt scattering evident in Fig 1 is per- 
haps surprising given the chiral nature of the pion and its interactions at low 




energy. Explaining this success will go a long way towards explaining the un- 
warranted success of the CQM in the light quark sector [28]. One way to do 
this is to construct a model of strong QCD which incorporates the physics of 
chiral symmetry breaking at a microscopic level. In this way one may compute 
TTTT scattering with composite particles while observing the dictates of chiral 
symmetry. It is possible to construct such a model by assuming a nontrivial 
QCD vacuum (typically a BCS-type vacuum) and building states on this vac- 
uum with the random phase approximation to the full Bethe Salpeter equation. 
This has been done in Coulomb gauge QCD [29, 30], with similar calculations 
in the Schwinger-Dyson approach [31] 

4 Conclusions 

Studying and understanding hadronic interactions is vital to hadronic physics 
and is an important part of nuclear and electroweak physics. It is probable that 
a microscopic description of hadronic interactions is necessary if one wishes 
to understand these phenomena in extreme conditions or in poorly known 
channels. The constituent quark model provides an excellent starting point for 
developing the understanding required to construct a reliable model of strong 
QCD. In the meantime, it also serves as an excellent phenomenological guide to 
the interpretation of scattering experiments. The development of continuum 
field theoretic models capable of describing hadrons and their interactions is in 
its infancy - we look forward to their maturation and application to reaction 
processes. 
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tal. A large portion of the work described here was performed with Ted Barnes 
and Adam Szczepaniak. The author is grateful to Paul Geiger for discussions 
on coupled channels. This work was supported by the US DOE under contracts 
DE-FG02-00ER41135 and DE-AC05-84ER40150. 
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